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Abstract 

The doubly differential cross sections for the 4 He(7r+, 7r~)4p reaction were calculated using 
both a two-nucleon sequential single charge exchange model and an intranuclear cascade code. 
Final state interactions between the two final protons which were the initial neutrons were 
included in both methods. At incident pion energies of 240 and 270 MeV the low-energy peak 
observed experimentally in the energy spectrum of the final pions can be understood only if 
the contribution of pion production is included. The calculated cross sections are compared 
with data. 
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1 Introduction 


The pion double-charge exchange (DCX) reaction is one of the rare forms of nuclear reactions 
which involve a minimum of two nucleons and provides a tool to study two nucleon correlations 
within a target nucleus. A number of studies have been made of exclusive DCX and interesting 
results have been obtained in regard to correlations[lj] and the propagation of pions through 
the nucleus in its ground state |3| 

In this paper we study inclusive DCX on a well understood nucleus where the principal 
physics interest lies in understanding and describing the mechanism. 

It is essential to be able to describe nuclear reactions on the hadronic level in order to 
interpret possible deviations from the expected behavior as “new physics” phenomena. Several 
techniques have been used to describe inclusive DCX. Becker and Schmit ||] used the sequential 
single charge mechanism in their calculations and their results showed general agreement with 
the energy dependence of the total cross section but predicted cross sections larger than the 
experiment data. 

Germond and Wilkin J5[] suggested a different mechanism for DCX reaction, the interaction 
of the incident pion with the exchange current meson or pions in the cloud surrounding each 
nucleon. Their calculations did not include the Pauli effect or Fermi motion of the nucleons. 
The production mechanism introduced by Jeanneret [|b] gave good agreement with DCX cross 
section at energies above the P 33 region and may be the main mechanism there. 

Gibbs et aZ.Q, reported another attempt to calculate the total and doubly differential cross 
section for the DCX reaction by using the sequential single charge exchange mechanism. They 
included a treatment of anti-symmetric wave functions. Both Pauli blocking and a crude form 
of final state interaction (FSI) were included in their work but interactions with the “spectator” 
protons were ignored as in the previous cases. 

In conjunction with dibaryon searches, Grater et al. [|j] calculated the total cross section 
for 4 He. In their model FSI was taken into account by using a Watson-Migdal Model in an 
eikonal approximation. 

Following the early total DCX cross section data, extensive experimental studies performed 
at Clinton P. Anderson Meson Physics facility reported high precision results for doubly dif¬ 
ferential cross sections. The experiments covered a wide range of elements (A=4, 208) over a 
range of incident energies (120 - 270) MeV. These studies showed a two peak structure in the 
energy spectra at forward angles at the higher energies in the low atomic number elements 
such as 4 He and 3 He |], 10], while that structure disappears in the heavier elements 16 O, 40 Ca 

ii- 

A number of calculations attempting to explain this two-peaked structure have followed 
a suggestion by M. Thies. He and Van Loon [IX] did calculations based on the fact that the 
angular distribution of the pion nucleon cross section is dominated by the p-wave and hence 
has a forward-backward peaked angular distribution in the 7 tN center of mass. In this case, 
pions coming from the DCX reaction in the forward direction would have made two small 
angle or two large angle scatterings. Two small angle scatterings would generate only a small 
energy loss while two large angle scatterings would bring about a large energy loss due to the 
recoil of the nucleons. This effect surely exists but there have been difficulties with obtaining 
a quantitative representation of the data using only this model. 

Another early attempt to explain the double peak structure in 4 He was by Wood||TT|| relying 


2 




on the mechanism proposed by Thies. Kinney || modified the Thies model by including several 
effects of the nuclear medium. His results gave reasonable quantitative agreement with the 
shape of the spectrum for 4 He for some angles, but did not follow the cross section as a function 
of angle. Recent calculations carried out by Kulkarni ||l4|| , based on a relativized version of the 
code of Ref. (7|] , found qualitative agreement with Kinney’s data. 

Part of the difficulty with the forward-backward argument is that the scatterings take place 
over a 3 dimensional sphere, which means that there is factor of sin# arising from the solid 
angle of the first charge exchange which mitigates the forward-backward angular distribution. 
Also, at 240 and 270 MeV, where the effect is the strongest, not only does the pion-nucleon 
angular distribution become more forward peaked in the center of mass as energy increases 
but it will also be thrown forward in the laboratory so that the backward (corresponding to 
low energy) peak is much diminished. We will return to this point in the discussion. 


2 Two-nucleon Model 

We first consider a model in which the two “spectator” protons do not play an active role in 
the scattering but do influence the result through the energy that they can carry off. 


2.1 Calculation of Matrix element 

We assume a sequential mechanism to describe the DCX reaction in 4 He(7r + ,7r _ )4p. The 


operator for the pion DCX amplitude in the plane wave limit [15] can be written as 


K D (k,k , ,ri,r 2 ) = 


1 

27T 2 


dqe 


—ik'-r 2 


/2(q, k' 


iq-(r 2 -ri) 

'—2—p-/ i( k, q) e 8kri 

q - K o 


where the pion-nucleon single charge exchange amplitude is represented as 
/(k, q) = [A 0 + Aik • q + i\ f a 2 ■ (k x q )]v(q)v(k). 

Here 2 1 2 is the pion propagator and v(q) is the off-shell form factor taken here to be, 

q — fc 0 


v{q) = 


A 2 + k 2 
A 2 + q 2 


( 1 ) 


( 2 ) 


(3) 


The quantity A (assumed here to be 4 frn 4 ) is related to the range of the pion-nucleon 


interaction and the A’s are obtained from 7T-N phase shifts| 1C ]. ri and r 2 denote the position 
vectors of the members of the struck nucleon pair. We first consider the double scattering to 
occur from two neutrons leading to a definite final state of 4 protons, as if it were an exclusive 
reaction. The density of these final states is then included to calculate the inclusive differential 
cross section. The matrix element for the reaction may be written as 


M = J dr 1 dr 2 dr 3 dr^ p (ri,r 2 ,r 3 ,r i ) < S’|K D |0 > ipHe(r i,r 2 ,r 3 ,r 4 ) 


(4) 


where if} 4 P ( r i, r 2 , r 3 , r,*) is the final state wave function for 4 protons. We assume that the four 
protons can be described by a product of plane wave states 
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(5) 


4 

V’4p( r i) r 2 , r 3 , r 4 ) = Y[e ipiri 
1=1 

'ipHei r) is the wave function of the initial state of 4 He which is taken here to have a Gaussian 
form: 

^He(ri,r 2 ,r 3 ,r 4 ) = Ne~ a ^ +S " +S $+*%, ( 6 ) 

where s, = r* — R, R = (iq + r 2 + r 3 + r 4 )/4, N 2 = (—) a -\/8, |0 > is the initial spin state 
for two neutrons in 4 He and < S'! is the final spin state of the struck pair. We assume a box 
normalization with unit volume. 

For the spin independent piece of the pion-nucleon amplitude 


where 


M 0 


g(P3,P ^ / dr 9 (r)e i ^~ p ) r e- ar2 
(27T 2 )2 J 


(7) 


with a 



1 (P3+P4) 2 (P3+P4) 

g(P3, P4) = 8 4 e 8a e 4a 


r p* qr 

g(r) = / dq-3- 2 f (q.) k 0/( k j h) 

J — q 

The differential cross section may be written as 


G?<r(k, k ; ) 

dQdk' 


\M\ 2 dp 3 dp 4 dQ p — 


( 8 ) 

(9) 


( 10 ) 


The integral over the angles of the final relative pp momenta is done by standard numerical 
methods while the 6 dimensional integral over p 3 and p 4 is performed by Monte Carlo by 
sampling the last factor in Eq. |§|. 

For the more general case of a fully spin-dependent amplitude, we must consider amplitudes 


Z 7 ' cr 2 i’ (T z 2 

r cr z i,(Tz 2 


=< a' zl ,a' z 2 \F D \a zl ,<j Z 2 > 


( 11 ) 


where Fd is the double scattering operator. The calculation of these amplitude is given is 
Appendix A. 


2.2 Final State Interaction 

The FSI is expected to play an important role in the low-energy region of the final nucleons 
(when the pion loses a small amount of energy) by increasing the value of matrix element. 
In this study, we introduce the effect of final state interaction among nucleons by taking into 
account of the interaction between the two protons in the pair which was produced from the 
two charge exchanges on the neutron pair. We calculate the correction coming from FSI by 
replacing the plane wave function for the active pair (the factor e _ip r in Eq. |?]) by <^(p) which 
contains an interaction in the s-wave. The modified matrix element including the FSI is 
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Mpsi 


( 12 ) 


where 


ff(P3,P4) 

(2vr 2 )l 


drg(r) 




<KP,r) 


</>(p,r) = 4irJ2i £ MPir)Y e m (r)Y e m *(p). (13) 

£,m 

Since we assume that the final pp interaction is important only in the relative s-wave 

Mfsi = !l |) '~ | P . : [ drg(r)e' lLJ i~ J - [0 o (p, r ) - jo(pr)\e~ ar2 + M 0 . (14) 

(27T 2 ) 2 J 

We approximate cf>o(p,r) as 


sin(pr + 5 q) — sin(<Jo)e anr 
pr 

with a n = 2ra T . 5q is the s-wave phase shift which is calculated by using the pp version of the 
(modified) Malfliet-Tjon potential^, [I^| . 

The strength of the FSI can be expected to have a strong correlation with the outgoing 
energy of the pion and produce the greatest effect at high pion energies. One can see this 
qualitatively by considering a DCX reaction on two isolated neutrons at rest. Since the total 
momentum of the recoiling pair of protons is equal to the change in the momentum of the pion 
we can write the conservation of energy as 

uj — <jJ = \J 4 (p 2 + m 2 ) + (k — k') 2 — 2m. (16) 

Thus, for a fixed angle, the momentum, p, is a monotonic function of the final pion energy 
and the larger the final pion energy, the smaller is the relative momentum of the two final 
nucleons and the stronger the final state enhancement. For this reason we can expect the 
inclusion of FSI to increase the cross section at high outgoing pion energies with little effect 
at low pion energies. Figure [l] shows the result of the two-nucleon model with (solid line) and 
without (dashed line) final state interactions. It is seen that the strongest effect occurs at high 
final pion energy where a significant increase in cross section is observed. 




3 Intranuclear Cascade Code 

To get more a realistic picture of the reaction dynamics, we need to include the protons in the 
initial state in our calculation which we do by using an Intranuclear cascade (INC) code in 
which the pion moves through the nucleus with probabilities of different kinds of interaction 
until it escapes or is absorbed. For an overview see Appendix B. 

Because the nucleons interact with classical potentials in this model, they can be deeply 
bound since there is no quantum constraint. To correct for this effect, these bound nucleons 
are freed and the binding energy removed from other particles. In the present case, bound 
nucleons were liberated with their final energy in their center of mass being a fixed number (20 
MeV) equal to about what they would have from their Fermi motion. This value was varied 
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and little effect was seen. The energy needed to free the particles was taken from the pion 
which is probably the right particle in many, but not all, cases. This feature is perhaps the 
most unrealistic aspect of the current INC code. This defect could perhaps be corrected using 
the method of Wilets et al. [19] but such a correction is beyond the scope of the present work. 


3.1 Final State Interactions 

The INC technique calculates the DCX cross section by following the classical motion of the 
pion through the nucleus. While there are interactions among all of the nucleons through 
potentials, as described in Appendix B, they are not of the quantum mechanical nature that 
was included in the two-nucleon model when FSI was treated. This QM correction is included 
in the INC by the use of a weighting factor for each event. When a DCX event takes place, all 
of the coordinates of the pion and the nucleons are known so that the matrix elements M$ and 
Mpsi can be evaluated. Since the event includes no QM FSI, one can correct the probability 
of its occurrence by including a weight which is equal to the ratio of the squares of these two 
matrix elements. 


3.2 Pion Production 

The aim is to model pion production, given a knowledge of pion production from a nucleon in 
free space. The total cross section for pion production from a free nucleon can be computed 
from the inelasticity in the phase shifts Jl^] since, for the moderate energies treated here, the 
(7r, 2n) reaction is the sole significant contributor to the inelasticity. 

Once it has been decided that a pion production event would take place in free space (based 
on the ratio of the pion production cross section to the total cross section), the distribution 
of the three-particle final state (two pions and one nucleon) in energy and angle must be 
modeled. To compute these functions the current relevant data are inadequate, so we have 
used two-step models. One possible model is that in which a A and a pion are produced with 
the A decaying to produce the second pion. This mechanism may indeed play a role. A second 
mechanism believed to be important just above threshold, would produce a “a” meson and 
then have it decay into two pions. The early data favors this mechanism, as does the fact that 
the production of an isospin zero pion pair is dominant. 

Here we use the “a” mechanism only in the sense that the pion-nucleon system is converted 
into a two-pion plus nucleon system (assumed isotropic in the original 7 tN center of mass) and 
then the two pion subsystem is allowed to decay (again isotropically). To model the appor¬ 
tionment of energy between the nucleon and the two-pion system we assume a distribution of 
the two-pion invariant mass proportional to Q(Q m ax ~ Q ) where Q is the momentum of one 
of the pions in their center of mass. Once these variables have been chosen, the momentum of 
the remaining nucleon is known. A Pauli blocking test (requiring that the nucleon must not 
have less energy in the final state than in the initial state) is applied, resulting in a significant 
fraction of the possible production events being disallowed. 

The total pion production cross section calculated at 240 MeV is 116 fib. Since the pro¬ 
duction cross section from a single free nucleon is about 600 fib, considerably less than the 
maximum possible cross section is obtained. This reduction is very likely largely due to the 
two protons which “shield” the neutrons from which production takes place. 
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3.3 Pauli Blocking 

The effect of Pauli Blocking becomes important at low energy. If there is little energy loss by 
the pion the energy of the 4 protons will be small. In this case there will be strong cancellation 
among in the Pauli terms as was discussed by Ref. [^]. Hence, we expect a reduction in the 
cross section due the Pauli principle. The effect of Paul blocking is taken into account in this 
study by using a weighting factor which consists of an incoherent summation of the squares of 
the expressions in Eqs. 0 and [18| (taken from Ref ]7j ). 


A] = 


1 


[24(1234) + 24(2134) - 4(1423) - 4(1324) - 4(24321) - 4(2341) 


- 4(3214) - 4(4213) - 4(3124) - 4(4123) + 24(3412) + 24(4321)] (17) 


A 2 = —^[4(1324)—4(1423)+4(2431) —4(2341)—4(3214)+4(4213)4(3124)—4(4123)] (18) 

v8 

where 


A(abcd) = e *(p- r i+Pb' r 2+Pc- r 3+p d - r 4) ( 19 ) 

The momenta in the first two positions correspond to those associated with the two nucleons 
on which the charge exchange occurs. Tests including this effect showed that it was very small 
at the energies considered here and this weight was not normally included in the calculation. 


4 Discussion 

We performed calculations for doubly differential cross sections by two methods, the first being 
described in Section 2 with the results shown in Fig. |l]. Our results at small angles clearly 
show the effect of the FSI which is around 30% at the peak at 25 degrees. 

It is clear from this figure that the two-nucleon model overestimates the absolute cross 
section in the region of the high-energy peak. The discrepancy could be due to the ignoring 
of scattering from the two “spectator” protons in our calculation. 

Figure ^ shows a comparison of the two models under the same conditions (i.e. only the 
two neutrons are included in the scattering in the INC). The same basic physics is present in 
the two models (one classical and one quantum) which are seen to give very similar results at 
240 MeV. 

Figure |3| shows the effect of including the two initial protons at 270 MeV. It is seen that 
the cross section is considerably diminished. The two protons are very effective in shielding 
the neutrons from the incident pion. 

Figure || shows the results with and without the FSI. It is seen that its inclusion by using a 
calculated weight in the INC is very similar to the direct calculation in the two-nucleon model 
(see Fig. [I]), the principal effect being at forward angles and high-energy outgoing pions. The 
shape of the high energy peak is influenced by the effect of the FSI. 
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These effects modify the cross section by up to a factor of 3. Figure j| shows the progression 
from the case with only 2 nucleons to the addition of the protons (without true absorption) 
to the addition of true absorption. 

Figure ^ shows the effect of FSI at 180 MeV. Our calculation overestimates the cross section 
but less than the other calculations i 0 • That could be due to the fact that 180 MeV is 
very close to the P 33 resonance peak, so that we expect a great deal of absorption to occur. It 
is also possible that the INC does not give good results at the resonance peak because there 
is a larger degree of coherence, not taken into account by the classical INC. 

While the full calculations (solid curves in Figs. [| and |5|) give good results for the high- 
energy peak, the low-energy peak is underestimated by about a factor of two at forward angles. 
Figures |?| and ^ show the effect of inclusion of pion production. Pion production complements 
the DCX channel at pion incident energies of 240 and 270 MeV, and influences the shape and 
value of the spectra in the region of the low energy peak. Vicente et aZ. [2C ] found that the pion 
production enhances the strength of the DCX cross section by 15% at 240 MeV and 30 % at 
270 MeV in 16 0 ( 7 r + , tt~)X. This same mechanism apparently explains the missing strength 
at the low energy peak at 270 MeV and 240 MeV in the present case. 

A measurement of pion production from 3 He was made in Ref. ]l(|. In this study negative 
pions were detected with a incoming ir + beam. Since there is no double charge exchange 
possible, the mesons must arise from pion production. Since there is only a single neutron 
in 3 He, whereas in 4 He there are two and since the blocking effects of the two protons should 
be about the same in the two cases, we expect the production in ^He to be about twice that in 
3 He. In Fig. [?] we show the data from Ref. © multiplied by a factor of two. The agreement 
with this data under this assumption is satisfactory. 

In regard to the possible forward-backward explanation of the two peaks, we note first 
that the effect is being observed in the laboratory. Even if the single charge exchange angular 
distribution remained forward-backward symmetric in the center of mass, it would become 
more forward peaked in the lab as the energy increases which would lead to the high energy 
peak becoming increasingly dominant. In fact the single charge exchange angular distribution 
becomes more forward peaked in the center of mass as well. These two effects lead one to 
expect a decreasing low-energy peak with rising beam energy since it is supposed to arise from 
the back scattering part of the angular distribution. The data show the opposite effect with 
the low-energy peak growing with respect to the high-energy one with increasing energy. This 
observation is naturally explained if pion production is important since its contribution rises 
with energy. 

Figure [9] shows the integrated total DCX cross section with the INC compared with data. 


5 Conclusion 

In this paper, we have presented calculations of the double charge exchange cross section for 
reaction 4 He (vr + , 7 r _ ) 2 p by using a sequential mechanism and an INC model. The considera¬ 
tion of FSI in our calculation enhances the high energy peak at small angles by 30 % leading to 
good agreement with experimental data. The missing cross section strength in the low energy 
peak in our calculation at energies 240 and 270 MeV is naturally explained by pion production. 
We note that there is a discrepancy between theory and experiment in regard to the shape 



of the 180 MeV data. While inclusive pion production gives only a very small contribution (one 
is only about 5 MeV above threshold) a significant contribution can come from the production 
leading to 3 He plus a proton. The threshold for this reaction is about 167 MeV. The calculation 
of this cross section is beyond the scope of the techniques used in the present paper. 

We wish to thank Mark Yuly for supplying the numerical values for the production cross 
section on 3 He from Ref. ©■ This work was supported by the National Science Foundation 
and the Jordan University of Science and Technology. 

A Spin Dependence 

To include spin in the calculation, it is preferable to rewrite the pion-nucleon scattering am¬ 
plitude in a slightly different form. 

By using the relations 


= q)rr'(q') 

v m,m' 


and 


8ni 


~ ~ / u/16 v > ^m.m — M-\/-m / * \~\rvrt / ~ 

’ q q = 37 ! ^ ^m , 1 (rAi (q) 


( 20 ) 


( 21 ) 


m,m' ,M 

we can then express the nN amplitude as 


f(q,q')= £ (^) 7 D»(i)q“q'”v«(q)v«(q') £ C”i”I,’ M t; M (i)Y”(q)Y” (q') (22) 

/3=0,1.2 v-a m.m'.M 

where a, (3 and 7 are related by the following table. 


p 

a 

7 

0 

0 

0 

1 

1 

0 

2 

1 

1 


The matrices to and t^ are defined by 


and 


+ o _ U) _ ^ . + i _ 1 _ 

1 1 — (7 — (J %, t -y — (J — 



[&x “t"” ^1 — ® — /— ^y) 


1 

72 v 


D 0 (i) = 47tA 0 (i); Di{i) 


,(i);0 2 (i) 


8m. 


The operator in spin space for the double scattering will then be given by 
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( 23 ) 


F D (k,k')= E ^ 1 (l)^ 2 (2)fc Q1 fc to2 y™ 1 (ki)F™ 2 (k 2 )x 

01,132 

V^/gl L, Q ll O!i,7l Ly a2,02,72 r 7i K 1 ) 1 ^ V Z 1 


xA m " m2 


iQ!! ,CX.2 


(ko,r) 


where 


A™?™ 2 (k 0 ,r) = / g“ i+Q2 i;7(q)y™ 2 (q) 


a 2 + k? 


o 


P*qr 


-r/q 


^a 2 + q 2 ) (fcg — q 2 ) 

We are now ready to calculate the matrix elements of this operator in spin space. Using 


we can write 


< a'\q M \a >= -V3C?~¥’ a ' 

2 2 


(24) 


<ai, 4 |F D (k,k')ki,(T 2 >= Fo^^Lk') 

= E (~ 1 ) Ml+M2D 0i( 1 ) D 02 (2)k“ 1 k , “ 2 F™ 1 (k)lh™ 2 (k')x 
01,02 

t ^ l mi,m' 1 ,Mi^ l m2,m l 2 ,M2 ( ^cn,—Mi,cr l 1/ ^ia2,-M<2,cr'2 A m i, m 2/i \ 

C 7 ai,ai, 7 i 070:2,0:2,72 01 ^ 1 > 7 i ^ 1 ^Kg,rj 


>,72, 


(25) 


B The Intranuclear Cascade 

The present code was originally developed to treat moderate-energy antiproton annihilation in 

^6|, |27|, I 2 HI • The annihilation 


nuclei and has been applied to that end several times [23, 24, [2^, £6], |27 
of an antiproton leads to pions (or at least it is so treated by the code) and so the history of 
pions in the energy range below and of the order of 1 GeV are essential to the calculation of 
energy deposition. It has had considerable success in predicting the rapidity distributions of 


strange particles produced in antiproton reactions[25]. The question of pion absorption and 


comparison with data has been addressed [29] and the code has been used for the comparison 
with inclusive data in work by Zumbro et al.\ 3C[. 

The initial bombarding projectile is started with appropriate initial momentum toward a 
circle which is large enough to contain the projected density of the target nucleus. A fraction 
of the beam particles (usually about one half) pass without interacting and cross sections are 
computed as the fraction of the reactions of interest which occur multiplied by the area of this 
disk. 


B.l Treatment of Nucleons 

The heart of the code is the construction of a model nuclear system with as many characteristics 
of a “real” nucleus as possible. The code exists in several versions according to the reaction 
being investigated and the detail required. 
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B.1.1 Version I: The Woods-Saxon Potential 

In this original model the nucleon-nucleon interaction is thought of as being made up of two 
parts: 1) a long range interaction (one pion exchange) and 2) a short-range part (heavy meson 
exchange, quark interaction, etc.). 

The average of the long-range part is represented by a potential with a Woods-Saxon shape. 
The parameters of the well are chosen so that the particle distribution that results from motion 
in this well represents the measured proton densities in nuclei. 

Each nucleon has a designated binding energy specified in the input data. For a given 
binding energy there is a point at the edge of the well where the potential is equal to the 
binding energy and hence the kinetic energy is zero. For each particle a limiting radius, R x , 
is chosen as that value of the radius and a set of position coordinates are chosen uniformly 
distributed within a sphere of radius R x . Once the position for a given nucleon is established, 
its potential energy can be calculated and, using the specified binding energy, the kinetic energy 
and hence the magnitude of the momentum, can be obtained. The directions of the momenta 
are chosen randomly in an isotropic manner. Since each particle is moving in a conservative 
potential well it will maintain the same total (binding) energy throughout its motion unless 
disturbed by an outside agent. 

The short-range part of the nucleon-nucleon interaction is represented by a scattering cross 
section. 


B.1.2 Version II: Inter-nucleon Potentials with Fixed Well 


In this version the potential well is generated by the superposition of semi-realistic nucleon- 


nucleon potentials due to the nucleons themselves, those by Malfliet and Tjon [17, 18]. 


If Vi is the potential seen by the i th nucleon at position r* then 


V i = JH^'(l r * _r jl) 


(26) 


where the subscript ij on V is to label whether the interaction is np or pp (nn). For the 


technique of creating the initial nucleus see section B.3. 

If we wish to guarantee that the nucleus remains in the ground state (i.e. with each 
nucleon keeping its assigned binding energy) we can fix the wells to those generated by the 
originally chosen nucleons with each nucleon moving in its own fixed well, even though the 
other nucleons which created this potential have moved to other positions as the time evolution 
progresses. Since the movement of non-struck nucleons is small, the actual wells should not 
be very different from these fixed wells. The reason for wanting to keep the fixed wells is to 
maintain the quantum condition that the nucleus remain intact with the same “quantized” 
states unless acted upon by an external force. As before, the nucleons can collide with a 
probability computed from a cross section entered as data. Note that this is necessary since 
the nucleons cannot exchange energy except by this mechanism. A nucleon may scatter from 
another nucleon’s initial position but it simply bounces off of the potential and the nucleon 
“struck” knows nothing of the interaction. 

One inconvenient point which is especially bothersome for light nuclei is that a single 
nucleon can be left bound in the phantom potential well. For example for 4 He if, in 7r + 
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absorption, three protons are removed, the remaining neutron is still bound in the well (if it 
is not specifically struck by the incident pion or one of the leaving nucleons) and is left with a 
negative energy when it should have a positive energy distribution (presumably related to its 
Fermi momentum). 

It is this version which was used for the inclusive measurements by Zumbro et al. 


B.1.3 Version III: Inter-nucleon Potentials with a Changing Well 

In this case the calculation of the potential wells follow the nucleons as they move. In some 
sense this is the most realistic of all the models. It has problems, however. As in any classical 
model the particles can evaporate from the system. If the nucleus falls apart rapidly the 
density at which the reaction takes place may not be correct, even though the initial nucleus 
is originally thrown with an acceptable density. To minimize this problem the nucleons do not 
begin to move until there is a first interaction of the projectile with one of them. Note that 
a typical nuclear reaction takes place on the order of a few fnr/c. Thus, while the nucleus is 
decreasing in density as a function of time, the problem is not very serious. 

B.2 Interaction of the Projectile with the Nucleons 

The next step is to allow the projectile to propagate through the nucleus and interact with the 
nucleons. Only pions will be treated in the present discussion. Kaons, ry’s and other mesons 
are also propagated by similar techniques. Immediately after moving the nucleons one time 
step the meson multiple scattering subroutine is entered. This subroutine moves the pion (or 
pions) one time step, taking into account the possibility that it can elastically scatter, charge 
exchange, be absorbed on a nucleon pair or produce one or more additional pions. 

For each pion the following steps are carried out. 

The distance from the pion to each nucleon is computed and compared with the value of 
R derived from the equation 


30|. 


&max — ItR (27) 

where a max is the maximum p ion-nucleon total cross section (that at the 33 resonance). If the 
distance is greater than R the index of the loop passes on to the next nucleon. The total cross 
sections used for comparison are 7r° cross sections, or in other words, the average of 7r + and 
7 t” cross sections. A weight is carried for the event which is the ratio of the true total cross 
section to the 7r° total. 

If the distance is less than R then the possibility that an interaction might take place is 
pursued. From the momenta of the pion and the nucleon the effective laboratory momentum 
for the pion is computed, i.e. the value the pion would have (for the same center of mass 
energy) if the nucleon were at rest. This is done because c.m. cross sections are tabulated as 
a function of Pi a b- The distance from the pion to the nucleon is then compared with a length, 
R, obtained from crriPiab ) = ^P 2 where ot is tabulated in a vector indexed by P[ a b- If the 
distance is greater than this number then there is no interaction with that nucleon. 

When a case is found that the distance to the nucleon is less than R then a weight for the 
event is computed from the ratio of cross sections and one of a series of branchings is selected 
according to conservation laws and the appropriate probability. 
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First, the possibility of pion absorption is explored. In order for a pion to be absorbed 
it must have interacted previously with at least one nucleon and absorption must be possible 
(charge conservation). If those two conditions are met, then a random number is compared 
with an input parameter, P a b s , to determine if an absorption will take place. 

If the pion is absorbed then the total energy of the pion is converted into momentum added 
to the current and previous nucleon in such a way that the sum of the additional momenta 
are zero (they are equal and opposite) and along the direction of the trajectory of the pion 
traveling between the two nucleons, i.e. parallel to the current pion momentum. If is the 
current value of the pion momentum then the current nucleon will receive 

Dl Pi_ 

P + — = «Ptt 

and the previous nucleon 

-P + ^=/?p ; 

To conserve energy we have (non-relativistically) 

Pi + aPvr) 2 , (P2 + /3Ptt) 2 _ P\ , P-2 , 

- ~T - — - I - \ OJ 

2m 2m 2m 2 m 

Where 1 (say) is associated with the current nucleon and 2 with the previous one. Solving 
for a and f3 the momentum P can be found. The direction of relative to the incident beam 
direction (in the case of only two nucleons) approximately follows a p-wave distribution from 
the previous elastic scattering. The direction of this large momentum dominates the angular 
distribution and tends to give a reasonably good representation of nd absorption. 

If pion absorption does not take place then the possibilities of pion production (if energy 
conservation permits), charge exchange (if charge conservation permits) or elastic scattering 
are chosen according to the ratio of the relevant cross sections to the total. 

If pion production does not take place then charge exchange or elastic scattering occur. In 
this case a table of Legendre polynomial coefficients are used (from a data table) to select an 
angular distribution for the final two particles. If the scattering or charge exchange leads to a 
nucleon whose energy is not greater than its initial value by the amount of the Pauli blocking 
energy the scattering is not allowed. 

B.3 Technique for creating a Nucleus with correlations 

The procedure starts from the point of view of the shell model. The full Hamiltonian can be 
written as 



(28) 

(29) 


+ rj) = E (31) 

i j^i 

This equation can be regarded as a classical or operator relationship. 

In the spirit of the shell model it is natural to write the total energy as the sum of single¬ 
particle energies E = Ej and make the association 

Ti + \ X! ~ r j) = E i ( 32 ) 

j 
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for each particle. Certainly this is not necessary but it is sufficient to satisfy the original 
equation and it allows us to associate an energy with each particle. These equations may 
be regarded as operator equations for the Schrodinger equation or simply classical energy 
conservation equations. A significant difference is that in the quantum case the kinetic energy 
symbols mean derivative operators and in the classical case they are non-negative quantities. 
In the quantum case they lead to negative multipliers of the wave function only in the case of 
tunneling. Since tunneling is a small probability we make the approximation that the kinetic 
energies should always be positive. In any event, we are forced to require that the kinetic 
energies be non-negative for the construction of the nucleus since it is to be used in a classical 
simulation. 

The resulting algorithm can be described as follows: 

1) All of the particles are assigned a value of binding energy in the shell model sense. As 
an example, for 12 C, there might be 4 particles with a binding energy suitable for the s-shell 
and 8 particles with a binding energy in the p-shell. The sum of these binding energies should 
be equal to the total binding energy of the nucleus. 

2) Each shell is assigned a sphere with some defined radius. This radius parameter allows 
a fine tuning of the radius for each shell to be obtained. If this radius is taken too small the 
nucleus will be compressed. If it is taken too large time will be wasted in attempting to create 
impossible configurations. In general it should be taken of the order of the effective square 
well radius of a given shell. 

3) A nucleon-nucleon potential is assigned to each particle pair. In general this is defined 
between particle types (nn and pp or np). 

4) The actual creation of the nucleus is begun with the throwing of a random uniform 
distribution of particles within the spheres defined above. 

5) The calculation of each kinetic energy is made using the positions of the nucleons to 
calculate the potential energy and subtracting it from the individual binding energies. If any 
kinetic energy is negative then the complete re-throw of the nucleus is started again with step 
4. If all kinetic energies are non-negative then the proposed nucleus is possibly viable. Note 
that the particles cannot be very close or the repulsive potential will make some kinetic energy 
negative nor can they be too far apart for the same reason. 

6) Now the magnitudes of the momenta are calculated from the kinetic energies and the 
directions are thrown in a random uniform manner. This process will lead to a nucleus with 
a non-zero total momentum. 

7) The total momentum of the nucleus is calculated by vectorially adding the individual 
momenta. If the magnitude exceeds a set limit, the angles of the momentum vectors are thrown 
again in an attempt to make the sum less than the limit returning to step 6. Since it is possible 
that no combination of angles will lead to a total momentum less than the prescribed limit, 
after some fixed number of tries (typically 300) the entire nucleus is re-thrown by returning to 
step 4. 

8) 1/A of the total momentum is subtracted from each nucleon thus putting the nucleus 
at rest. This procedure does not leave the binding energy unchanged but if the limit on the 
total momentum is small the error on the binding energy is also small. For example, it was 
found that if the total momentum is restricted in step 7 to be less than 40 MeV/c. then the 
maximum error in the binding energy was observed to be less than 0.2 MeV (out of 28 MeV) 
for 4 He. 
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Figure 2: Comparison of the INC with only two nucleons (neutrons) active (histogram) with the 
two-nucleon model (solid line). 
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Figure 3: The results of INC model at incident energy 270 MeV with only the two neutrons (dashed 
histogram) and four nucleons (solid line histogram). 
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Figure 4: The results of doubly differential cross section for the reaction 4 He(7r + , 7r - )4p at incident 
pion energy of 240 MeV. The calculation of the solid histogram includes FSI and the dashed line 
histogram does not (absorption is included in both). 
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Figure 5: The results of the INC at incident pion energy of 240 MeV. The solid histogram represents 
the results including FSI and absorption factor, the short dashed histogram has no absorption. The 
long dashed histogram represents the results of taking two neutrons only. 
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Figure 6: The same as Fig. (| but at incident pion energy 180 MeV. 
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Figure 7: The results of doubly differential cross section at 240 MeV. The solid histogram shows our 
results with FSI, absorption with pion production, while the dashed line shows the result without 
pion production. The doted line shows the pion production cross section alone. The data given in 

m 


the solid squares are the production cross sections from Ref, 


multiplied by two. 
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Figure 8: The same as Fig. 8 but at an incident energy of 270 MeV. 









Figure 9: Total DCX cross section for 1 He(7r + , 7T“)4p. The solid line shows the results including 
FSI and absorption factors, while the dashed line is without FSI. The experimental data were taken 
from m (solid circles), 0 (open triangles), |HJ (open squares) and [Q (solid squares). 
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